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Abstract 

A Bose-Einstein condensate of cold atoms is a superfluid and thus responds to rotation of its 
container by the nucleation of quantized vortices. If the trapping potential is sufficiently strong, 
there is no theoretical limit to the rotation frequency one can impose to the fluid, and several phase 
transitions characterized by the number and distribution of vortices occur when it is increased from 
to oo. In this note we focus on a regime of very large rotation velocity where vortices disappear 
from the bulk of the fluid, gathering in a central hole of low matter density induced by the centrifugal 
force. 



1 Introduction 
1.1 The model 

The main theoretical framework for the description of rotating BECs is the so-called Gross-Pitaevskii 
(GP) theory. The behavior of the fully Bose-condensed gas is described by a single macroscopic wave 
function ^ and the ground state of the condensate is determined by the minimization of the GP energy 
functional. Assuming a strong confinement along the rotation axis so that the condensate is quasi-2D, 
the GP functional can be expressed, in the rotating frame, as 

^p'iysl*] - |i |(iV + A) + {v,,M ' l^'iotr^) 1*1' + ^} d'-- (1) 

Here A — iliotesAr with ^l^ot > the rotation velocity, 63 the unit vector in the xs-direction, r = (xi, X2), 
r = |r|, Vtrap the trap potential and l/e^ with e > the GP coupling constant. We work in units such 
that h and the particle mass are both 1 and the wave-function is normalized as / j^'p = 1. 

In ([T]) we have introduced a 'vector potential' A in order to separate the contribution of Coriolis and 
centrifugal forces due to the transformation to the rotating frame. This emphasizes the importance of 
the effective potential 

K=ff(r) - ytrap(r) - ir!?,,r2, (2) 
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and makes it immediately clear that the physical behavior of the gas will be strongly affected by the type 
of trap that is being used. In this note we consider a potential of the form 

^trap(r) = kr^ + \niy (3) 

with k > and s > 2 [anharmonic traps). Indeed in the purely harmonic (quadratic) case k — 0, 
the rotation frequency would be limited by the frequency fiosc of the trap and rapid rotation regimes 
forbidden: Any larger rotation speed would make Vcs unbounded from below, and thus the atoms would 
fly apart under the action of the centrifugal force. If A: > on the other hand, the effective potential 
is bounded below for any firot, and there is thus no theoretical obstruction to the exploration of regimes 
of very large rotation frequencies. 

It is also worth mentioning here that among anharmonic traps (s > 2) one has to distinguish between 
'soft potentials' (s < oo) discussed here and 'hard potentials' or 'flat traps' considered in [H [71 [11] 
(formally given by s = oo). Indeed, the physics of the transition to a giant vortex state is different in 
these two cases and the third critical speed is considerably smaller in the latter (after an appropriate 
scaling, see below). 

1.2 Phase transitions and critical speeds 

Let us briefly review what is expected to happen when the rotation velocity imposed to a trapped Bose 
gas is increased. We only discuss the case of anharmonic traps {k > 0) and refer to the literature (see 
e.g. [31 [3] for reviews) for the specificities of the purely harmonic case. 

• As long as the rotation speed is smaller than a first critical speed ftc^ , the fluid is vortex-free. 

• When flrot is increased past fici , vortices start to appear and form regular patterns. When many 
vortices are nucleated they are distributed according to a triangular lattice. This behavior lasts 
until a second critical speed flc2 is attained. 

• At fiiot ~ 1 the centrifugal force becomes so important that it dips a hole of strongly depleted 
matter density at the center of the trap. A vortex lattice survives in the annular bulk of the fluid 
until a third critical speed is reached. 

• When r^i-ot crosses the third critical speed ftcs, all the vortices retreat in the central hole of low 
matter density. The bulk of the fluid is then vortex-free but a phase circulation remains around 
the hole, indicating that it acts as a multiply quantized giant vortex. 

Experiments using anharmonic traps we are aware of have so far been limited to rotation frequencies 
r^iot < ^^C2 ■ In a clear density dip is observed at the center of the trap but experimental limitations 
have prevented the creation of the giant vortex. The theoretical problem is nevertheless of interest and 
has been intensely investigated P [TUl [TTl [T ^ [T ^ [Ti l [T5 ] . 

In this note we focus on the third phase transition, when the state of the condensate changes from 
a "vortex-lattice-plus-hole" to a "giant vortex". We explain the main physical insights leading to a 
mathematically rigorous estimate of the critical speed we have recently obtained. A more complete 
discussion and the details of the proof can be found in [51 [B] , together with a more thorough discussion 
of the literature. 

2 Emergence of the giant vortex state 
2.1 Scaling of the GP functional 

In the rest of this note we will consider the strongly interacting (Thomas- Fermi) regime e — >■ and focus 
on angular velocities firot ^ 1 close to the third critical speed. 



2 



In order to make the effect of the quadratic part in the potential Vtrap apparent even in the limit 
J7rot — i> cxD, we keep the ratio fJeff/i^rot fixed by setting O^g = j^'^^t some < 7 < 1 fixed, where 
fleS = (rifot ~ ^osc)^^^) so that the effective potential becomes Vcs{r) = kr'^ — iri^^r^. 

Moreover we scale variables in such a way that the bulk of the mass of the condensate does not expand 
as firot — > 00: Most of the mass is indeed concentrated close to the minimum of the potential 'I4ff(r) at 
f — Rm with i?ni = ['^Is/ isk)]^^'''^~^\ which wanders to infinity as iljot ~> 00. It is thus convenient to 
define the scaling in such a way that the global minimum of the rescaled potential is independent of the 
angular velocity. This can be done by introducing the following variables 

x = i?^ir, ij{x) = Rn,^{r), = iJ^^^ot (4) 

which yields f^ysi^"] = Rm^^'^^H'] with the scaled energy functional 

S'^^m - + ^^e^W + -fn^VixM^ + £-2|^|4| dx (5) 

and the potential (with minimum at x = 1) 

V{x) = ^x' - ^x^. (6) 

The potential and interaction terms in ([S]) become comparable when ^ which is the order of 
the second critical speed fic2 (see [HIS]). In the sequel we take » e^^ (and will in fact specialize to 
n > e"*) and investigate the features of the transition associated with the third critical speed. 



2.2 The matter density profile 

As a first step we look for a good approximation to the density profile of the fluid. To this end we define 
the functional 

fgv[5] - [l^gf + l^'ix - x-'fa' + 1^' (^-/ - l^') 9' + ^-^/} dx. (7) 

This is nothing but ([S|) restricted to wave functions of the forn{3 

7/;(x) = .g(x)exp(il7i9) (8) 

with a real valued function g, normalized such that J g^ — 1. In the giant vortex ansatz ([5]) all the 
vorticity is concentrated in a single vortex at the origin with huge winding number fl. We denote ^gv 
and Egv the ground state and ground state energy of £gv respectively. 

A noticeable fact is that when ft — iloe^^ with Hq — 0(1), (?gv progressively changes from a Thomas- 
Fermi to a gaussian profile. This is most conveniently understood by a change of variables. Since the 
minimizer of ([7]) is radial we consider only radial functions and define y = il^^^{x — 1), and f{y) = 
Q~^/^g{x). Equation (O then turns into 

£:-[/] ^ + laVf + ^o'^'f] dy (9) 

up to an unimportant constant term. We have set = 4 + 7(,s — 2), Taylor-expanded the potential and 
neglected terms beyond quadratic (see O Section III.B]). Without the last term the minimizer is the 
gaussian 

/o.c(2;) = 7r-i/4ai/4exp{-iay2} (10) 

^Strictly speaking we should use the integer part of f2 in the phase factor, but since £7 is very large in our regime, it 
makes little difference to assume that it is an integer. 
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which win give a good approximation (after scaling) to the giant vortex profile g^v when is large, i.e., 
O ^ e~^. Notice that it is quite different, in particular because of its long tails, from the TF profile 
(with /xtf a chemical potential fixed by normalization and [ • ]+ the positive part) 

fMy) = f^y''^[MTF-i«V]+, (11) 

which is a better approximation of g^^ for ^ e^"'. 

An important step of our approach is an energy decoupling: Defining, for any normalized wave function 
tp a function v by writing ip — ffgvW exp(ir2i9) and using the variational equation satisfied by t^gv, we get 
£:GP[^] ^ E'^'' +£[v] with 

= fflv {^|V«r - B ■ J(t;) + £-^(1 - |t.|')'| dx. (12) 

We denote by J(t;) — ^{v\/v* — w*Vv) the superfluid current associated to v and B = 0(a; — x~^)e^ the 
vector potential corrected by the contribution of the giant vortex. 

Our analysis of the phase transition leading to the giant vortex state is based on the minimization 
of the reduced functional (fT2)) under the constraint / ^gvlt'P — 1- Thanks to the energy decoupling 
this is equivalent to the original variational problem: A function u minimizes (jl2p if and only if ip — 
5gvU exp(if2'i?) minimizes the scaled GP functional ([S]). Note that ggv does not vanish, except at the 
origin, and thus all the potential vortices of ip are accounted for by the reduced wave function u. 

2.3 Disappearance of vortices 

Now that we have a good approximation to the matter density profile, we can discuss the nucleation of 
vortices in the bulk of the fluid, that is in the region where ggy is significantly large. We first discuss 
the size that vortex candidates would have. Using an ansatz v that vanishes inside the bulk in a region 
(vortex core) of size t and minimizing ([7]) with respect to t gives an estimate of the size of vortex cores 
: t ^ s^/^fl~^/^. From the considerations in the above section, one can guess that the thickness of the 
annulus where the mass resides is of order fi^^/^. When fl > e^^ we have t > il^^/^ which means that 
vortices become too large to completely fit in the bulk of the fluid, whereas they are relatively small when 
^l -C This alone is of course no proof that a phase transition occurs in this regime but it turns out 
that one can, when fl <i; e^*, rigorously estimate the energetic contribution of vortices and prove that 
many are nucleated and uniformly distributed in the bulk. This result is obtained by making use of the 
powerful "vortex balls methods" (see, e.g., [IJII3 for reviews) that apply in this regime. 

The regime > e^^ on the other hand poses a new challenge. Since vortex candidates are comparable 
in size to the bulk of the fluid, the intuition behind the vortex balls methods fails and a different approach 
has to be used. In fact, one can rule out density depletions without any precise estimate of the vortices' 
energy. The argument goes as foUowfH: a trial state v — 1 indicates that if u minimizes (dH), then 
£[u] < 0. If on the other hand we can prove that £[u] > with a sufficiently small error, it is clear that 
it is not favorable to nucleate vortices, since the trial state w = 1 is then almost optimal. 

To estimate the energy it is convenient to integrate the second term of p2p by parts : Since V • (7gyB = 
0, we can write (?gvB — V^F with the dual gradient = {~dx2, dx^) and a suitable potential function 
F. By Stokes' formula we then hav(j^ 

- / gl^B-Jiu)^ [ FV^-3{u). (13) 

^With significant simplifications, see [SJE] ^^'^ ^^^^ details. 

•^In full rigor, the integration should be restricted to a bounded region and boundary terms included in the analysis. 
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Now, it is clear that the only potentially negative term in (IT^ is ([T^ . It is thus the only one that could 
favor the nucleation of vortices. Note also that the quantity • 3{u) appearing in the r.h.s. of ([T^ is 
nothing but the curl of the superfluid current. It can thus (in analogy with fluid mechanics) be thought 
of as a vorticity measure. A trivial estimate yields | • J(u)| < |Vwp and we thus have 

/ ^fflvlVup-/ glB-3{u)> f fLl^-\F\)\\7u\\ (14) 

The potential F is an explicit function of ^gv and, using for the latter the gaussian ansatz described 
in the preceding section, one can see that 2|i^| < g^^. When JIq is sufficiently large for the gaussian 
approximation to become accurate one can then infer that (I14p is > 0, which indicates that nucleating 
vortices can not decrease the energy. This is the desired lower bound to the energy E\u\ and actually 
gives more information: Since we have not used the third term in (jl2p so far, we deduce an upper bound 
to this term. A careful analysis then reveals that energy considerations prevent u from vanishing in the 
bulk of the condensate. 

As mentioned before the behavior in the case of hard trapping potentials (s — cxd) is quite different: the 
squeezing of the condensate against the hard walls of the trap due to the centrifugal force is much more 
pronounced is this case and vortices therefore energetically more costly. This explains why the third 
critical speed is considerably larger for soft trapping potentials (fica ^ £~^) than for hard potentials 
(17c3-^e-2|loge|-i). 

3 Conclusions 

Wc have sketched the energetic considerations leading to our proof that a phase transition from a vortex 
lattice to a giant vortex state occurs in BECs rotated at a very large rotation speed. Our main finding 
is that in the rescaled variables introduced in Section l2.ll the transition happens when f2 3> i.e., 
^rot ^ e-4(s-2)/(s+2) jj-^ -j-j^g physical variables. It would be interesting to know a precise estimate of the 
critical speed for the phase transition, but since the potential vortex cores are comparable in size to the 
bulk of the fluid, it is not even clear that a sharp transition, rather than a smooth crossover, occurs. 
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